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С помощью разработанной численной методики рассматриваются способы понижения значений параметров пространственно- 
го трёхмерного напряжённо-деформированного состояния ведущего устройства сборки в зоне соединения его секторов и ме- 
таемого стержня. Реализуется подход, основанный на ужесточении этой наиболее нагруженной части секторов с помощью вкла- 
дышей, выполненных из прочного материала. Полученные результаты позволяют ставить и рассматривать вопрос о некоем по- 
добии рационального проектирования геометрии вкладышей секторов сборки. Задача упругопластического деформирования 
решается с помощью вариационно-разностного метода. Физические соотношения принимаются согласно теории малых упруго- 
пластических деформаций. Г еометрические соотношения берутся в виде уравнений Коши. Физически нелинейная задача реша- 
ется методом переменных параметров упругости. 


1. Введение 

Для анализа напряжённо-деформированного 
состояния (НДС) метаемых элементов и вспомога- 
тельных, т.н. ведущих устройств (ВУ), с помощью 
которых осуществляется разгон сборки при высо- 
коскоростном метании, в [1-4] представлены чи- 
сленные методики расчёта двумерного осесимме- 
тричного и трёхмерного пространственного термо- 
упругопластического деформирования конструк- 
тивных элементов сборки с использованием квази- 
статического подхода, а также результаты числен- 
ного счёта и анализа их НДС. 

Решение проводится вариационно-разностным 
методом, реализующим экстремальный вариаци- 
онный принцип Лагранжа для функционала пол- 
ной потенциальной энергии системы «тело-на- 
грузка» методом конечных разностей. В общем 
случае функционал энергии является неквадратич- 
ным. Для анализа упругопластического поведения 
конструкций используется деформационная тео- 
рия пластичности А.А. Ильюшина, которая хоро- 
шо описывает процесс деформирования при моно- 
тонном нагружении и широко применяется в рас- 
чётах упругопластических тел. Решение физически 
нелинейной задачи сводится к решению последо- 
вательности квадратичных задач с уточняемыми в 
каждом приближении параметрами нелинейности. 
С помощью этих параметров осуществляется дви- 
жение вдоль диаграммы нелинейной зависимости 
интенсивности напряжений от интенсивности де- 
формаций <т~е„ которая аппроксимируется лома- 
ной двухзвенной линией. Для минимизации ква- 
дратичной функции многих переменных использу- 
ется условие экстремума сеточного аналога функ- 
ционала энергии, которое сводит проблему к реше- 
нию системы линейных алгебраических уравнений 
большого порядка относительно искомых переме- 
щений в узлах конечно-разностной сетки. Матри- 
ца системы уравнений имеет ленточную структуру, 
симметрична и положительно определена, что 
очень важно при численной реализации. Решение 
получающейся системы уравнений для двумерных 
задач проводится прямым методом Гаусса, для 


трехмерных задач используется итерационный ме- 
тод верхней релаксации с выбором оптимального 
коэффициента релаксации. 

Работы [3, 4] посвящены численному модели- 
рованию пространственного трёхмерного термо- 
упругопластического НДС выполненного в форме 
катушки ВУ при метании тяжёлого недеформируе- 
мого стержня массой 2 кг (рис. 1). Причиной про- 
странственной постановки задачи является кон- 
структивная трёхмерность ВУ, состоящего из трёх 
одинаковых дюралюминиевых секторов РРNНСМ 
(рис. 2), между гранями которых находится предот- 
вращающий прорыв газов герметик. К левому тор- 
цу ВУ и левой части стержня (рис. 1) приложено 
давление газа Р, посредством которого сборка раз- 
гоняется в трубе метательной установки (МУ). За- 
дача решалась в декартовой системе координат х, у, 
I. Геометрия секторов ВУ задавалась следующими 
параметрами: I - длина, Д, - внешний радиус; К 2 - 
внутренний радиус (радиус стержня); К 3 - радиус 
цилиндрической части. Описанная модельная за- 
дача реализовывалась для следующих значений 
геометрических параметров: Д,=5,0 см, К 2 = 1 ,0 см, 
Г=15,0 см. Ввиду равноправности секторов при на- 
гружении, а также симметрии параметров НДС от- 
носительно плоскости Охі, при численном анализе 
(как и представлении полученных результатов) 
рассматривалась половина ММ НО сектора с углом 
<р=60°. Соединение секторов ВУ со стержнем осу- 
ществлялось по части контактной поверхности 
«стержень-сектор ВУ» (т.н. гребёнка) с радиусом К 2 
и значениями осевых координат 2 ГР Н и 2 ГР В , ограни- 
чивающими соответственно слева и справа область 
гребёнки. Таким образом, в рассмотренной поста- 
новке задачи гребёнка представляет собой зону 
приложения нагрузки со стороны метаемого 
стержня к исследуемым секторам ВУ. 

При метании сборки для значений давления 
Д=-250 МПа, модуля упрочнения материала 
Е Х =Е / 5 и радиуса цилиндрической поверхности 
Д 3 =3,5 см в верхней части центрально расположен- 
ной по длине сектора гребёнки длиною І ГР =4,0 см 
(2’ ГР Н =5,5 см, 2 ГР В =9,5 см) возникают растягиваю- 
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щие, а в нижней её части - значительно большие по 
величине сжимающие осевые напряжения. Ампли- 
туда изменения этих напряжений вдоль гребёнки 
составляет Дсг г =1180 МПа. Наибольшие значения 
величин интенсивности деформаций и напряжений 
находятся на этом же участке внутренней поверхно- 
сти г=К г секторов: (е / ) тах =0,04 и (сг,) тах =665 МПа со- 
ответственно, и достигаются в нижней части гре- 
бёнки вблизи края грани для значения угла ср=54,4°. 
Для значения модуля упрочнения Е Х =Е / 25 интен- 
сивность деформаций (е,) шх =0, 18, что является не- 
допустимо большой величиной для данной кон- 
струкции. Таким образом, проведённый анализ по- 
лученных напряжений и деформаций во всём объе- 
ме секторов как по длине, так и по их толщине по- 
зволил выявить зону гребёнки «...как наиболее на- 
груженный локальный участок секторов, где отме- 
чается значительная неравномерность распределе- 
ния напряжений и деформаций, способная приве- 
сти к разрушению сборки. Причиной такого ра- 
спределения параметров напряжённо-деформиро- 
ванного состояния на гребёнке является нагрузка от 
метаемого стержня» [3]. Использование в качестве 
материала секторов сплава В95 с более высокими 
прочностными характеристиками (предел теку- 
чести оу=589 МПа, деформация начала текучести 
е=0,0085) привело к значительному снижению ве- 
личин параметров НДС ВУ. Так, величина макси- 
мальной интенсивности деформаций по сравнению 
с дюралюминиевыми секторами уменьшилась при- 
мерно в три раза: (е,) шах =0,014. Всё это позволяет 
увеличить давление газа, а следовательно, - и ско- 
рость метания сборки. Тем не менее, уже для значе- 
ния давления />=400 МПа наибольшая величина 
интенсивности деформаций в секторах (е,) шах =0,153. 




Рис. 2. Г еометрические параметры секторов ведущего 
устройства. ~ один из трёх секторов ; МЫНС 

- рассчитываемая половина сектора 


Для устранения неравномерности распределе- 
ния и понижения уровня напряжений и деформа- 
ций на гребёнке и во всём ВУ вместе с применени- 
ем более прочных материалов также были предло- 
жены и численно реализованы достаточно очевид- 
ные технические решения [4]: 1) увеличение пло- 
щади контактной поверхности «стержень-сектор 
ВУ» за счёт увеличения длины гребёнки / ГР ; 2) уве- 
личение толщины цилиндрической части секто- 
ров, определяемой радиусом К } . Так, для сплава 
В95, давления />=400 МПа, длины гребёнки 
/ ГР =7,0 см и Д 3 =3,5 см величина интенсивности де- 
формаций (е,) шах =0,018. Для / ГР =4,0 см величина 
(е,) гаах =0,097 для радиуса /? 3 = 4,0 см и (е,) тах =0,04 для 
радиуса 7? 3 =4,5 см. Все эти меры позволили снизить 
величины параметров НДС на гребёнке секторов 
до значений, позволяющих гарантировать не толь- 
ко штатное функционирование сборки, но и безо- 
пасность процесса метания. Тем не менее, следует 
заметить, что увеличение массы сборки (на 7 и 16 % 
соответственно) по сравнению с первоначальной 
массой 4,066 кг (для /? 3 =3,5 см) за счёт изменения 
толщины цилиндрической части секторов приво- 
дит к уменьшению скорости метания сборки. Из 
представленных результатов видно, что вариант 
увеличения длины гребёнки является наиболее 
перспективным при решении задачи уменьшения 
параметров НДС в секторах. Проведённые расчеты 
для случая / ГР =/ и /?,=3,5 см для сплава В95 позво- 
лили без изменения массы сборки увеличить вели- 
чину давления газа до значения />=550 МПа при 
наибольшем значении интенсивности деформаций 
в секторах (е,) тах =0,05. 

В случае фиксированных геометрических раз- 
меров секторов ВУ, места расположения гребёнки в 
сборке и её длины, значительно меньшей, чем дли- 
на секторов, изложенные подходы и способы 
устранения неравномерности распределения пара- 
метров НДС в зоне гребёнки и уменьшения значе- 
ний этих параметров, как это видно из предста- 
вленных результатов, не позволяют находить прие- 
млемое решение. Поэтому необходимо искать 
иные пути решения задачи: например, использова- 
ние для ВУ более прочных и лёгких материалов; 
новые конструктивные решения для секторов и 
т. д., имея при этом в виду нежелательным увеличе- 
ние массы метаемой сборки. Ниже предлагается 
один из таких подходов. 

2. Расчёт НДС секторов ВУ с вкладышами 

Для уменьшения значений параметров НДС ВУ 
при метании сборки предлагается наиболее нагру- 
женную часть секторов в зоне гребёнки ужесточить 
с помощью вкладышей, выполненных из прочного 
материала. На рис. 3 с соблюдением пропорций для 
описанных выше геометрических размеров пред- 
ставлена половина сектора ВУ с вкладышем 
(область, расположенная в центральной части сек- 
тора и закрашенная тёмным цветом). Вкладыши 
одинаковы по форме и размеру, а каждый из них в 


131 




Известия Томского политехнического университета. 2006. Т. 309. № 6 


отдельности представляет собой полученную двумя 
осевыми сечениями (для значения угла между ни- 
ми <р= 120°) треть толстостенной трубы длиной Г вкл 
и радиусами К 2 и /? вкл (см. рис. 2) внешней и вну- 
тренней поверхности соответственно. Длина вкла- 
дышей определяется значениями осевой коорди- 
наты и 7 ВКЛ В . Гребёнка, размеры и положение 
которой задаются осевыми координатами 
2 ГР Н =5,5 см, 2 ГР В =9,5 см, имеет длину І ГР =4,0 см и 
расположена симметрично относительно значения 
осевой координаты %=Ь/ 2, т. е. в средней по длине 
части секторов. На рис. 3 гребёнка обозначена 
областью, закрашенной по сравнению с вклады- 
шем более светлым тоном. 



Рис. 3. Рассчитываемая половина сектора АР5МЫНС с вкла- 
дышем 

Задача решается при следующих статических и 
геометрических граничных условиях для половины 
сектора [3]: 

• торец г=0: 

ст = ст = 0 , а = Р: 

• торец і=Ь: 

ст = ст = ст = 0: 

Х2 у 2 2 ’ 

• контактная поверхность Г/ѴЯ.Ѵ «стержень-сек- 
тор ВУ», кроме гребёнки, (жёсткая стенка со 
скольжением): 

и = ѵ = 0, 

"У ", =° (Я, =0); 

• гребёнка (жёсткая стенка без скольжения): 

и = ѵ = 0, 

<7 „ И.Ѵ + ѴугПу = Р ст (<т„ = Р С1 ); 


поверхность АМИР (у=0) (рис. 2, 3): 
ѵ = °, ст^ = а г = 0; 

контактные поверхности «труба МУ-сектор ВУ» 
А В, ЕМ (жёсткая стенка со скольжением): 
и = ѵ = 0, 

а :Л+°у-Пу = ° (О‘ в =0); 

свободные конические поверхности ВС, БЕ (ра- 
венство нулю проекций на оси координат дей- 
ствующего на наклонной площадке полного на- 
пряжения): 

+с ѵг +а М\ = °> 

+<Т , П у +<Т уЛ =0: 

а =Рх +а * п у +а -л =°; 


• свободная поверхность СВ, задаваемая ради- 
усом Ку — 

+ °\,«ѵ + 2 Сду «, ф = 0 (а,. = 0), 

(а у - о х ) п х п у + а ху (п х -Пу) = 0 (ст, „ = 0), 

+я,«, =° (я, =о); 

• свободная поверхность грани В НО, для которой 
(р=60° (рис. 2, 3) - 

СТ Х +о-Х _ 2гт "■ ". = 0 (<% = °)> 

К ” а -ѵ ) "ѵ »ѵ + % А - К ) = 0 (ст,. , = 0), 

- сг _-х« > = 0 (<Г-„ = °)> 


где сг ѵ , о у, а ѵ ст,, а уѵ ст,, а ух , ст,, ст г , - компоненты 
тензора напряжений; и, ѵ, к — компоненты вектора 
перемещений в декартовой системе координат; п х , 
п у , п у - направляющие косинусы; ст,, ст ѵ , ст,, ст„ ? , ст ж , 
ст,, ст,,,, ст.,,, ст,, — компоненты тензора напряжений в 
цилиндрической системе координат г, <р, 21, Р С т _ 
приложенная на гребёнке нагрузка от стержня. На 
осевые перемещения \ѵ ограничения не накладыва- 
ются. Для наглядности рядом с некоторыми стати- 
ческими граничными условиями в скобках содер- 
жатся записи этих условий в цилиндрической систе- 
ме координат. На контактных поверхностях вклады- 
ша задаются условия непрерывного контакта. 

Результаты получены на конечно-разностной 
сетке (/х/х&)=( 15x17x66) при реализации системы 
примерно 50500 линейных уравнений итерацион- 
ным методом верхней релаксации. Здесь і - коли- 
чество узлов сетки по толщине (по радиусу), у - по 
окружности, к - по длине секторов. Материал сек- 
торов - сплав В95. Вкладыши выполнены из высо- 
копрочной хромоникельмолибденованадиевой 
стали 38ХНЗМФА, имеющей высокие предел теку- 
чести ст=1080 МПа и величину деформации начала 
текучести е =0,0051. Рассматривается вариант на- 
гружения сборки давлением Р= 400 МПа. 

Как упоминалось выше, при отсутствии вклады- 
шей наибольшая величина интенсивности дефор- 
мации в секторах (е,) тах =0,153, а величина 
(ст,) т х=889 МПа. Эти значения достигаются в одной 
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и той же ячейке на внешних боковых гранях 5НС ря- 
дом с внутренней поверхностью г=Я 2 в средней по 
длине части сектора. Максимальная амплитуда осе- 
вых напряжений в ВУ Д<х=1715 МПа имеет место 
также на гребёнке: в верхней её части для 2=9,4 см 
<х=695 МПа и <х=— 1020 МПа в нижней части для 
2=5,4 см. В силу характера действующих внешних 
нагрузок осевые напряжения изначально являются 
определяющими в НДС метаемой сборки. Практи- 
чески во всем объеме сектора ВУ эти напряжения 
сжимающие. Исключением является небольшая по 
размерам область в верхней части гребёнки на кон- 
тактной поверхности «стержень-сектор ВУ», огра- 
ниченная по длине сборки значениями 7,5 и 10,0 см 
осевой координаты и 1,5 см по радиусу, где осевые 
напряжения являются растягивающими. Около 
1,5 % объёма ВУ деформируется пластически. Ниже 
гребёнки для значения осевой координаты 2*4,6 см 
материал секторов достаточно близок к выходу в 
область пластического деформирования по всей их 
толщине (от внутренней поверхности радиусом Я 2 
до внешней поверхности радиусом Л,). При увели- 
чении величины давления в секторах ВУ возможно 
образование т.н. пластического шарнира. В напра- 
влении осевой координаты размеры этой зоны нез- 
начительны, т. к. они представлены одним рядом 
ячеек конечно-разностной сетки. 

Учитывая полученные результаты для наиболь- 
ших значений интенсивности напряжений и де- 
формаций, место положения вкладышей в секто- 
рах было выбрано симметричным относительно 
значения осевой координаты 2=С/2. На рис. 4 
представлены зависимости максимальных значе- 
ний интенсивности деформаций (е,) шах от толщины 
к и длины Х вкл вкладыша. Толщина вкладыша опре- 
деляется соотношением: к=Я вш -Я 2 . Кривые 1, 2, 3 
построены для случая симметричного расположе- 
ния вкладыша в секторе ВУ для трёх вариантов зна- 
чений длины І вкл , равных 1,5; 3,5 и 5,0 см соответ- 
ственно. 

Как видно из представленных на рис. 4 результа- 
тов, наличие вкладыша приводит к уменьшению 
значений параметров НДС, в данном случае - вели- 
чины интенсивности деформаций. При этом следу- 
ет отметить комплексный характер влияния геоме- 
трических параметров вкладыша на этот процесс: 
чем больше длина Х вкл вкладыша, тем значительнее 
уменьшение интенсивности деформаций с увеличе- 
нием толщины к по сравнению с вариантом без 
вкладыша. При выборе геометрии вкладышей 
необходимо иметь в виду, что увеличение их длины 
и толщины приводит к увеличению массы секторов 
ВУ и, как следствие, - к уменьшению скорости ме- 
тания сборки. Поэтому выбор геометрии вкладыша 
следует делать по результатам всесторонней оценки 
её влияния на НДС сектора с учётом изменения 
массы сборки. В рассмотренных вариантах увеличе- 
ние массы сборки за счет наличия вкладыша для 
минимальной и максимальной толщины к соста- 
вляет 0,22 и 1 ,94 % для варианта 1; 0,53 и 4,50 % для 


варианта 2; 0,76 и 6,50 % для варианта 3 соответ- 
ственно. Наибольшие величины максимальных 
значений интенсивности деформаций (е,) шах имеют 
место на внутренней поверхности г=Я 2 секторов как 
в нижней части вкладыша, так и ниже вкладыша в 
зоне гребёнки. Величины (е,) та „ начиная, пример- 
но, со значения 0,017 и меньше, достигаются на 
внешней цилиндрической поверхности 

Анализ полученных результатов показал, что 
верхняя часть вкладышей менее нагружена по срав- 
нению с её нижней частью, как это было в зоне гре- 
бёнки секторов при отсутствии вкладыша. Поэтому 
был рассмотрен вариант несимметричного относи- 
тельно значения координаты 2=С/2 расположения 
вкладыша, для которого на рис. 4 построена кривая 
4. В этом варианте 2 ВКЛ Н =5,00 см и 7 ВКЛ В =8,75 см, а 
его длина Х вкл =3,75 см. Увеличение массы сборки 
для минимальной и максимальной толщины к со- 
ставило 0,56 и 4,90 % соответственно, что (для слу- 
чая максимальной толщины) на 1,6 % меньше, не- 
жели в варианте, представленным кривой 3, в кото- 
ром длина вкладыша Х вкл =5,0 см. Таким образом, с 
точки зрения наименьшего увеличения массы ме- 
таемой сборки вариант несимметричного располо- 
жения вкладыша является предпочтительным по 
сравнению с вариантом симметричного его распо- 
ложения относительно координаты 2=С/2. 



О 0,2 0,4 0,6 0,8 1,0 

Рис. 4. Зависимости максимальных значений интенсивности 
деформаций ( е,Д„ в секторе ВУ от геометрических 
параметров вкладышей 

Представленные на рис. 4 результаты позволя- 
ют выбирать геометрические параметры вкладыша 
для получения заранее заданного значения макси- 
мальной интенсивности деформаций (е,) гаах в секто- 
рах ВУ. Например, для получения в секторах НДС, 
интенсивность деформаций при котором не пре- 
вышала бы величины (е,) шах =0,033 (это значение на 
рис. 4 отмечено пунктирной горизонтальной лини- 
ей), можно воспользоваться следующими вариан- 
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тами геометрии для вкладыша: 1) к=к А ~ 0,125 см, 
^вкл = 5,0 см; 2) к=к^ 0,14 см, І вкл =3,75 см; 
3) к=к(~0,68 см, І вкл = 3,50 см. 

Что касается вероятности появления пластиче- 
ского шарнира, то она, судя по уменьшению объе- 
ма пластически деформируемого материала секто- 
ра, уменьшается. 

3. Выводы 

На основе разработанной численной методики 
рассмотрены способы уменьшения значений пара- 
метров пространственного трёхмерного напряжён- 
но-деформированного состояния ведущего устрой- 
ства (ВУ) сборки в зоне соединения секторов ВУ и 
метаемого стержня - т.н. гребёнки. Проведённые 
ранее исследования, основанные на использова- 
нии в качестве материала секторов лёгкого и проч- 
ного сплава В95 и достаточно очевидных конструк- 
тивных решений — увеличении площади контакт- 
ной поверхности «стержень-сектор ВУ» за счет из- 
менения длины гребёнки и толщины цилиндриче- 
ской части секторов, - позволили ценой значи- 
тельного повышения массы метаемой сборки уме- 
ньшить интенсивности деформаций. Отмечается, 
что эти подходы (особенно второй) для сборки дан- 
ной геометрии практически уже не позволяют по- 
вышать скорость метания за счет увеличения да- 
вления газа. 
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Предложенное и численно реализованное в ра- 
боте конструктивное решение для секторов сбор- 
ки, имеющее целью ужесточение наиболее нагру- 
женной в зоне гребёнки области с помощью вы- 
полненных из прочного материала вкладышей, да- 
ло положительный эффект. Полученные результа- 
ты позволяют ставить и рассматривать вопрос о не- 
коем подобии рационального проектирования гео- 
метрических параметров вкладыша секторов ВУ 
метаемой сборки. Для этого необходимо построить 
семейство кривых зависимости (е,) тах =ДХ вкл , А) для 
разных вариантов расположения вкладыша (и не 
только для вкладыша рассмотренной цилиндриче- 
ской формы) с целью получения номограммы, 
пользуясь которой, для заданного значения макси- 
мальной интенсивности деформаций (е,) таі можно 
определить геометрические параметры вкладыша. 

Относительно созданной методики следует от- 
метить, что она позволяет получать информацию 
обо всех параметрах пространственного напряжён- 
но-деформированного состояния с точностью до 
ячейки конечно-разностной сетки при проектиро- 
вании и оценке прочности не только рассмотрен- 
ного варианта сборки, но и других конструкций 
для широкого спектра внешних нагрузок. 
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